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Abstract 



We show that processes involving pions that remain very near their mass- 
shell can be reliably computed in the pionless effective field theory, with the 
pion integrated in as a heavy field. As an application, we compute the vr- 
deuteron scattering amplitude near threshold to next-to-leading order in the 
momentum expansion. This amplitude is formally dominated by an infrared 
logarithm of the form log (7/7717^), where 7 is the deuteron binding momentum, 
and rn-TT is the mass of the pion. The coefficient of this logarithm is determined 
by the S-wave pion-nucleon scattering lengths. 



I. INTRODUCTION 



For many processes involving pions scattering from multi-nucleon systems, the kinematics 
are such that the pions must be included as dynamical particles, as they can carry energy 
and momentum that is not small compared to their rest mass, and can be far off-shell. 
However, in the case of low-energy pion-deuteron (vrrf) scattering one has a situation where 
the pions can carry energy and momentum that puts them on or near their mass-shell, both 
on external legs and internally in some diagrams. Furthermore, given the small binding 
energy of the deuteron, the two nucleons in the deuteron are off their mass-shell by an 
amount that is far less than the pion mass. Therefore, we can construct an effective field 
theory (EFT) where non-relativistic nucleons interact via contact operators, and pions near 
their mass-shell are included as massless particles that generate a coulomb-like potential. 
The construction of this low-energy EFT is similar to the development of the heavy-quark 
effective theory (HQET) §g and non-relativistic QCD (NRQCD) ^-Q. 

The EFT relevant to momentum transfers much less than the pion mass, in which nu- 
cleons interact via contact operators, is known as the pionless EFT, or EFT(7f). The power- 
counting rules in EFT(7f) are now an established part of the nuclear physics literature and 
can be found in several review articles [Q, together with reviews of the substantial array of 
phenomena that have been studied in EFT(7f). A recently developed variant of EFT(7f) is 
particularly suited to high-order, precision calculations 0. In this EFT, known as dibaryon 
pionless EFT, or dEFT(7f), the ^Si NN effective range parameter, r^, is taken to be of order 
l/Q in the power counting, where Q is a small momentum This naturally leads to the 
use of dibaryon fields to account for nonperturbative enhancements 0. In dEFT(7f), opera- 
tors are ordered according to powers of the total center-of-mass energy [0. Unlike EFT(7f"), 
the higher-dimension operators involving external fields are not renormalized by the S-wave 
strong interactions, and therefore do not scale with inverse powers of the renormalization 
scale. Thus, naive dimensional analysis of these operators is sufficient to estimate their con- 
tribution to a given process. Due to these technical advantages, we will work with dEFT(;f) 
in this paper. 



II. INTEGRATING IN THE PIONS 

Given that both EFT(7f) and dEFT(7f) have been studied in detail, the 

remaining ingredient is the pion. In analogy with HQET and heavy-baryon chiral pertur- 
bation theory (HBxPT), for processes involving the pions very near their mass-shell it is 
convenient to perform field redefinitions to remove the classical trajectory of a pion and use 
nonrelativistic kinematics. In the limit that the nucleon mass is infinitely larger than the 
pion mass, and the deuteron binding energy is taken to be much smaller than any other 
scale in the problem, the pion kinetic energy in the deuteron rest frame, T^, can be taken to 
be a constant even in loop diagrams. Furthermore, the theory can be constructed in terms 
of the annihilation and creation operators for the pions. In loop diagrams the typical pion 
residual kinetic energy is of order while its momentum is of order 7, where 7 is the 

deuteron binding momentum. Thus, as is the case in NRQCD, it is possible and necessary 
to treat the residual energy operator in perturbation theory. Working in the isospin limit, 
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the leading-order (LO) Lagrange density describing non-interacting "heavy pions", Hh, is [| 



£ = tr 



(1) 



where is the kinetic energy of the pion and where is the pion mass. The subscript 
"/i" denotes a heavy field, and we have defined the meson matrices 
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(2) 



There is an implied summation over all possible pions with kinetic energy T.^, just as there 
is a summation over all possible four-velocities in HQET. The operator vr^ annihilates a pion 
while the operator vr^ creates a pion. The propagator for the pion has the form 



^T.(g°,q) 



(3) 



2m^T^ — |q 

which is independent of and reduces down to that of a "potential photon" in NRQED or 
a "potential gluon" in NRQCD, (i.e. — i/(|qp — ie)) in the limit that the pion has vanishing 
kinetic energy, — 0. 

We describe deuteron properties using dEFT(7f), in which the scattering length and 
effective range in the ^S'l-channel are assumed to be unnaturally large, and to scale like 1/Q 
in the power-counting The Lagrange density describing the dynamics of the nucleons 
and dibaryon, , in the '^5'i channel is 

V2 



ida + 
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N - tl 



V2 

ido + - A 
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t] N^P^N + h.c. 



(4) 



where y is the coupling between nucleons in the ^5*1 channel and the ^^i-dibaryon. The 
spin-isospin projector for the ^Si channel is 



Tr 



p«t pj 



(5) 



It is easy to show that this Lagrange density alone reproduces the NN scattering amplitude 
in the ^5'i-channel, with scattering length and effective range ra, when 



where /i is the renormalization scale. This identification is made after dressing the dibaryon 
propagator, as shown in Fig. 0. Higher order contributions to the scattering amplitude 
arising from the shape parameter and higher order terms in the effective range expansion 
are understood to be included perturbatively. 



^We could have chosen a different normalization for iTh such that = \/2m^7r/i, in which case 
the kinetic term would be 

v2 ^ 



£ = 7r 



+t 



TT^ + ... 



which is a more familiar form for the non-relativistic Lagrange density. The ellipses denote the 
kinetic-energy terms for the vr^ and vr'^. 
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FIG. 1. The dressed dibaryon propagator. The bare dibaryon propagator is dressed by nucleon 
bubbles to all orders. Each diagram counts as in the power- counting scheme. 




FIG. 2. Diagrams with interaction vertices that do not conserve pion number (left) appear 
as contact interactions (right) in the heavy-pion EFT. The single (double) dashed lines represent 
ordinary (heavy) pions and the solid black lines represent nucleons. The solid bar with a blob is a 
dressed dibaryon field. 



In order to have the incoming and outgoing nucleons near their mass-shell and also to 
have all the pions in a process near their mass-shell, there must be an equal number of 
pions entering and leaving each interaction vertex. As illustrated in Fig. 0, processes with a 



different number of pions entering than leaving [p!5|-[T7[| will be represented by local operators 
(on the scale of m^) in the EFT. One can assign a charge to the pions -not carried by the 
nucleons- which must be conserved in all processes in the EFT. Interactions between a heavy 
pion and a nucleon are described by a Lagrange density of the form 



C = co(T^) A^^A^ tr 



+ ci(T,) A^V'^A^ tr 



(7) 



As there are no loop contributions from the heavy pion fields ^, the coefficients Co,i(T^) 
represent the complete pion-nucleon (irN) scattering amplitude evaluated at the center of 
mass energy T^. Hence, these coefficients are related to the parameters in the effective range 
expansion for nN scattering via 



-47r(l+/?) 



1 



k cot (5^jv ~ iV^'ITT'nTn 



^For the bubble contributions to the vr/jA^ scattering amplitude, one can close the contour in the 
upper half of the complex energy plane to avoid enclosing a pole in the energy integration. Thus 
all bubble contributions vanish. This generalizes to all bubble diagrams involving heavy pions. 
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( +'^) cot 6^j^ — iy/2mj,Tj, ' ^ ^ 



where the effective range expansion for k cot 6^j^ has the form 



1 



A;cot5±v = + - ' (9) 

where are the scattering lengths in the s-wave isoscalar, (+), and isovector channels, (—), 
and Sttn are the phase shifts for vrA^ scattering. The ellipses denote the contributions from 
the effective range terms and higher. We have retained the quantity f3 = mj^/MN in these 
expressions as it arises from a trivial kinematic relation between the scattering amplitude 
and the Lagrange density when T-n- = 0. However, it is important to keep in mind that 
the there will be contributions from higher orders in the expansion that have factors of f3 
associated with them due to the explicit expansion. At threshold, these coefficients 

reduce down to 

Co(0)= 47r(l + /3) a+ , Ci(0) = -27r(l + /3) a' . (10) 

To determine the coupling between the heavy pions and the dibaryon field in the ^5*1- 
channel, we realize that in the limit that the pions only interact once with the nucleons, 
and using the fact that the deuteron is isoscalar, the only contribution to the ird amplitude 
is from the operator with coefficient Cq. Furthermore, this contribution must be absolutely 
normalized, as the pions couple to the baryon-number operator N^N. This gives rise to a 
uniquely determined interaction between the heavy pions and the dibaryon field, described 
by 

C = - 2co(T^) t] tr [MlMh] . (11) 

As the interactions induced by the heavy pions are subleading in the momentum expansion, 
it is anticipated that the modifications to eq. ([l^) occur only at subleading order. Indeed, 
at NLO in the EFT one expects a contribution from a counterterm whose coefficient is not 
dictated by baryon-number conservation. We write this contribution as 



mImA , (12) 



which scales like ~ Q in the dEFT(7f) power-counting. Here A ~ and /i is the renormal- 
ization scale. In general, the coefficient function ?7(T7r,/i) will have both real and imaginary 
parts, as bubble diagrams involving the pion vanish. The factors of M^v, A and enter as 
per the power-counting appropriate for dibaryon fields P. 



III. Tid SCATTERING 

At low-momentum, the LO amplitude for vrrf scattering receives contributions from the 
diagrams shown in Fig. ^ There is the one-loop diagram with a single insertion of the 
T^N scattering amplitudes in eq. (|^, and the local ttc? operator given in eq. (pi]). At LO, 
the operator that is inserted in the deuteron is a charge operator, and consequently it is 
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(a) 



(b) 



FIG. 3. Graphs contributing to ird scattering at LO. The dark ovals in diagrams (a) and (h) 
are interaction vertices taken from eq. ^ and from eq. (1J_), respectively. 



absolutely normalized at zero-momentum transfer, and the form factor is identical to that 
of the deuteron charge form factor |^ at LO. Therefore, at LO, the nd scattering amplitude 
is 

4 /Ikl 



|k|) = 2 co(r.) 



7^3 



tan 



(13) 



where 7 = \/M]\fB is the deuteron binding momentum, with B the deuteron binding energy. 
The momentum transfer to the deuteron is k. As expected, this form factor reduces to 
A^^f^ = 2 co{Tt,), proportional to the isoscalar charge of the deuteron, in the limit k — >• 0. 

At next-to-leading-order (NLO) in the EFT expansion, there is a contribution from a 
two-loop diagram (rescattering) involving a heavy-pion propagating between the nucleons, 
and from a local nd operator, as shown in Fig. El The amplitude at NLO is 




(a) 



(b) 



FIG. 4. Graphs contributing to n-d scattering at NLO. The logarithmically divergent two-loop 
graph arising from heavy pion exchange, diagram (a), is renormalized by the local ird operator given 
in eq. H, denoted by the solid square in diagram (b). 



7 



(1 - 7r3) 



( 



27r 



co(r^)'-8ci(T^) 



log 



4(27-zV2m^T^)% |k| 



+ 



4(27-^72^^) 



tan 



-1 



2 27-iV2m^T^ 



(14) 



where we have used the integrals given in the Appendix of Ref . [|T^ . The two-loop diagram 
is logarithmically divergent and is regulated in the IR by the deuteron binding energy for 
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incident pions with vanishing kinetic energy. The /i-dependence of the counterterm, "tiITt^, /i), 
must exactly compensate the /i-dependence of the loop diagram. We have used the relation 
TjiTT,,^) = r7(T^,/i)/(AM7v). Given that the scale where dEFT(7f) is expected to breakdown 
is A ~ m^, we choose to renormalize the theory at /i ~ m^r. 
In the Tjr, |k| — limit this amphtude reduces to 



ANLO) 



7 



7'^3) 



7^(0, /i) - - [co(0)2-8ci(0)' 

TT L 




(15) 



It is important to notice that there is a contribution that is non-analytic in the deuteron bind- 
ing momentum, that is logarithmically enhanced, ~ log(7/m7r), arising from the two- loop 
diagram. This logarithm was found previously in the work of Borasoy and GrieBhammer []TU[ . 

The scale dependence of the logarithm determines the naive- dimensional analysis (NDA) 
estimate for r/(0, m.„): 



~ - |co(0)^ 
7r 



^Ci(O)^ 



(16) 



Note that the logarithmic divergence in the two-loop diagram modifies the NDA analysis for 
the dibaryon operators of Ref . . It is straightforward to show that other contributions to 
the amplitude, both momentum-dependent and momentum- independent, are higher order 
in the power-counting. Therefore the sum of amplitudes given in eq. ([T3|) and eq. ([M]) 
constitute the complete amplitude up to NLO. 

While we have constructed the EFT appropriate for describing the very-low energy be- 
havior of 7id scattering, there is additional information that we have at our disposal from 
HBxPT. We know the chiral expansions of the coefficients Cq and Ci; at LO, Ci(0) ~ m^r 
while Co(0) ~ [^. Armed with this information we see that, in fact, A^^f ^ ~ Co(0) 



and W"^^^'^ ~ ci(0)^ are of the same order. This does not represent a breakdown of the 
low-energy EFT, but results from the LO amplitude being anomalously small. Thus, given 
the anomalous behavior of the LO amplitude, the leading contribution to low-energy nd 
scattering is 



A 



nd 



ALO) 



ANLO) 



+ 



(17) 



as given in eq. (0) and eq. (0). The ellipses denote terms that are higher order in the 
momentum expansion, suppressed by additional factors of 
To NLO, the nd scattering length is given by 



dnd 



(1 + /5/2) 



;i + /?) a 

27 
1 - 



1 + (a+)2-2(a-)2 log 



27 



^(0,^7r) 

IGvr 



(18) 



It is interesting to compare this formula with that given in Ref. [^. There, instead of the 
logarithm, one has the overlap of the coulomb-like pion propagator between deuteron wave 
functions. Of course, since logarithms are universal, the long-distance tail of any sensible 
deuteron wave function will generate the infrared logarithm. The choice of different deuteron 
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wave functions corresponds to a different choice for the counterterm, fj{0, m.^). In principle, 
one may be able to determine 7^(0, m^) uniquely using wave functions computed in a 
fully consistent pionful EFT [g|%6|. 

We can estimate a^^ by taking f/(0, m^) = 0, which corresponds to choosing the scale of 
the logarithm to be set by 171.,^. We use the Neuchatel-PSI-ETHZ (NPE) pionic hydrogen 
measurement |^ f\ of the isovector vrA^ scattering length, a~^p = (—0.0905 ± 0.0042) m~^. 
The isoscalar vrA^ scattering length, = (—0.0022 ± 0.0043) is consistent with 

zero and therefore we will neglect it in our estimate of a,rd- Using m,r = ''^71+ = 
139.57 MeV and Mn = (M„ + Mp)/2 = 938.92 MeV, we find a^d = -0.019 m-^ as compared 
to the NPE pionic deuterium measurement ||2^ of Re (a^^) = (—0.0261 ± 0.0005) m'^. 
Therefore, we see that the leading non-analytic piece provides the dominant contribution 
to TTO? scattering. However, as expected, the contribution from the counterterm fj is an im- 
portant contribution that cannot be discarded. It is worth pointing out that an imaginary 
part of a^f will be induced at higher orders in the expansion. This is consistent with the 



,-1 



28 



smallness of the measured value of Im (a'^f) = (—0.0063 ± 0.0007) 

The nd scattering amplitude can be expressed entirely in terms of the vrA^ and nd scat- 
tering lengths, and the difference fj{T.,r,m.,r) — 77(0, m^). At the order to which we are 
working, this latter difference will be completely imaginary, and will be determined from 
matching conditions. In particular, the scattering amplitude must be unitary (we are not 
allowing for energy transfer from the pion) and thus the optical theorem will determine 
fj{T.^,m.^) — 77(0, mjr), since fj is independent of the three-momentum transfer. Moreover, 
this difference will be independent of the renormalization scale at this order. Expanding out 
the scattering amplitudes in 



v^27?v7V/''^7r, we find 



AUT., |k|) = 47r(l + (3/2) a^, + ^^ii±^a 
.87r(a+)2v/2^^(l + /?)7r3 



+ 
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1 - 7r3 

4(27 - W2m^T^ 
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4 (27 - zV2m^T^) +|k| 



+ 2 
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(27- 





(19) 



This expression provides the formally leading energy and momentum-transfer dependence of 
the nd scattering amplitude. The formal EFT expansion parameters are 'y/m.jr and pTr/m-jr, 
where p.^ is the pion momentum. Since 'y/rriT^ ~ 1/3, we expect important higher order 
contributions, arising from the momentum dependence of the vrA^ scattering amplitude (and 



^We do not use the determinations of the scattering lengths in Ref. 
ment of Re {a-j^d) 



which include a measure- 
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its effective range expansion) and also from higher-order nd local operators involving two or 
more derivatives, which are suppressed by additional factors of 

In order to determine r7(T^, m.„) — fj^O, m^) we integrate the square of the LO amplitude 
to obtain the cross section for 7vd scattering, and then enforce unitarity order by order in 
the Cj's. It is convenient to define the function, S{x), where 



S{x) 



4 (tan ^ X 



logx 



TT 



- i4tan-^x ( Lz2(e*'*""''") - Li2{-e 
+ 2 ( Li3(e^2*^°"'^') - Li^i-e'^ ''''''' ^) 



i2 tan 



2 / -1 1, 

a; tan x log 

irs V 2 



1+x' 



X 

Y 



(20) 



which, by inspection, is real valued, and at x = has the value 5(0) = ( ^ ^^^^ ) ■ The 
functions Lin{z) are the polylogarithmic functions of order n. To the order we are working, 
unitarity implies 



f]{T^, m^) - 77(0, m^) = zStt 




7^3 



S 



, 27 , 



tan 



V 27 , 



(21) 



In this expression we have not included contributions of order (3 or higher as they are higher 
order in the expansion. Evaluating eq. ([T9|) numerically and using eqs. (pOD and (|2l|) we find 
that the amplitude at threshold is Xd(0, 0) = - 0.46 fm. While at = 20 MeV we find 

Ad(20,0) = ( -0.30 -^0.00013 ) fm , Ad(20, 150) = ( -0.27 - i0.038 ) fm . (22) 

Clearly the imaginary part of the amplitude is subleading in the EFT expansion. 



IV. CONCLUSION 

In this work we have shown how to include pions into the pionless effective field theory. 
There are kinematic regimes in processes involving more than one nucleon with pions in the 
external states, where there are infrared enhancements resulting from pions being near their 
mass-shell in intermediate states. We have constructed an effective field theory to describe 
these processes and have analyzed nd — >• nd. The scattering amplitude for this process is 
dominated by the leading non-analytic contribution of the form ~ log (7/m7r), which we can 
simply compute in this new effective field theory. 
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